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ABSTRACT 


A   fundamental,    non-linear,    ordinary  differential   equation   is 
derived    for    the    flexural  vibrations    of  a    non-uniform  beam,    including 
shear. 

This   equation   is    solved  with    the    specializing   assumptions    of  a 
constant    shear   coefficient   and    that    the   moment    of   inertia    is   directly 
proportional    to   the    cross    sectional   area. 

A   perturbation  method    is    used   and    the   eigen   frequencies   and  mode 
shapes    obtained   are    functions    of   a    perturbation   parameter   e.      These 
solutions   are   accurate    up   to   the    second    order   of   e. 
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SYMBOLS   AND  UNITS 


ENGLISH    LETTER   SYMBOLS 

A  Beam  cross   sectional  area,    1?" 

Lengthy   combinations    of   other   quantities    (see   equations    (14) 
and    (14a) 

Proportionality   constant,    I   =  C..A,    L 

-2 
Young's    modulus,    FL 

-2 

G       Shear  modulus,  FL 

4 
I       Moment  of  inertia,  L 

k  Shear  coefficient,  dimens ionless 

L  Beam  length,  L 

M  Bending  moment,  FL 

s  kAG(y'-\|r),  F 

t  Time,  T 

V  Shearing  force,  F 

x  Distance  along  beam,  L 

y  Beam  deflection  (positive  downward),  L 

GREEK  LETTER  SYMBOLS 
o*       LT  ,  radians 

©     TT^j  L 

2   -4 
p  Mass  density,  FT   L 

cr    kiq,  f"1 

\lr  Slope    of   beam  centerline   due    to  bending 

CO  Frequency   of  vibration,    T 


OTHER 

SYMBOLS 

/ 

A 

• 

_4_ 

Symbols  not  seen  here  may  be  found  in  Appendix  A. 
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INTRODUCTION 


The  Timoshenko  beam  equations  have  been  well  established  [5]*. 
Neglecting  rotary  inertia,  these  equations  become: 


(i) 

(2) 


The  left  hand  sides  of  both  equations  (1)  and  (2)  are  functions  of 
y  and  i|r .   Because  of  this  common  dependence,  these  equations  are 
referred  to  as  coupled.   Uncoupling  of  equations  (1)  and  (2)  involves 
their  reduction  to  one  equation  in  which  either  y  or  'if,    but  not  both, 
appear. 

It  has  been  shown  that  the  correction  for  shear  is  much  greater 
than  that  for  rotary  inertia  [10].   The  effect  of  shear  becomes 
increasingly  important  as  the  frequency  of  vibration  is  increased 
[9].   Therefore,  the  effects  of  shear  are  worth  investigating. 
Figure  (1)  is  a  free  body  diagram  of  a  beam  element.   Note  that 
the  quantity  (|bt-^)  in  Figure  (1)  is  the  slope  of  the  be; 
centerline  due  to  shear  deformation. 

Figure  1 


:am 
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*List  of  references  appears  on  page  53 


The    shear   term   in  equation    (1)    contains    the   shear   coefficient,   k. 
Timoshenko  defined    the    shear   coefficient   as    average    shear   stress    over 
maximum  shear   stress.      This    shear   coefficient,    then,    is   a    function  of 
the   cross    sectional   area    of   the   beam.      Timoshenko   found    it    to   be   con- 
stant   for  general   classes    of  geometrical   shapes.      However,    there    have 
been   other   studies    investigating    this    shear   coefficient   which   have 
disagreed   with   Timoshenko's   values.      In    1966,    G.    R.    Cowper    [4] 
derived   expressions    for    the    shear   coefficient   which   were    functions 
of    the   geometry,    Poisson's    ratio,    and    the    dimensions    of    the    cross 
section    in    some    instances.       R.    D.    Mindlin   and    H.    Deresiewicz    [9] 
maintain   that   Timoshenko's    shear   coefficient    is    a    function  also   of 
the    frequency   of  vibration.      The    following    table    contains    some 
representative   values    for    the    shear   coefficient.      These   values 
were    determined    by  different   methods    and   different    assumptions. 
Therefore    their   comparison   is   difficult. 

Table    1    (From  Cowper    [4]) 

k    for   rectangle  k    for   circle 

Timoshenko  0.667  0.750 

Mindlin  0.822  0.847 

Cowper        u   =  0  0.833  0.857 

u  -   .3  0.850  0.886 

u  ■   .5  0.870  0.900 

There    seems    to  be   a   conflict    of   opinions    in   the    literature    on   the 
shear   coefficient.      This    paper   develops    a   general  equation  with  k 
as    a    function   of   beam    length   but    all    solutions    treat   k   as    a    con- 
stant.     This    is    justified    by    the    fact    that    for    the    class    of   areas 
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investigated,  that  is,  where  the  moment  of  inertia  is  directly 
proportional  to  the  cross  sectional  area,  there  is  general  agreement 
that  k  is  a  constant. 

The  remainder  of  this  paper  will  be  devoted  to  the  uncoupling  of 
equations  (1)  and  (2)  and  the  subsequent  solving  of  the  uncoupled 
equation  for  its  eigen  frequencies  and  associated  mode  shapes.   This 
uncoupled  equation  will  be  general  with  respect  to  cross  sectional 
area,  moment  of  inertia,  and  shear  coefficient.   That  is,  these 
quantities  will  be  functions  of  the  distance  along  the  beam. 
Solutions  of  this  type  would  be  of  use  in  optimization  of  design. 
It  is  often  desirable  to  minimize  deflection,  bending  stress  or  modify 
frequency  of  vibration  in  beam  design.   By  allowing  the  cross 
sectional  area  of  the  beam  to  vary  with  x  it  may  be  possible  to 
reduce  deflection  and  bending  stress  or  to  modify  the  frequency  of 
vibration  through  appropriate  selection  of  an  area  function.   In 
this  thesis  the  area  function  will  be  defined  in  such  a  way  that 
only  one  dimension  of  the  cross  sectional  area  be  a  function  of  x. 
For  the  sake  of  comparison  to  existing  solutions,  this  variable 
dimension  must  vary  in  such  a  way  as  to  keep  the  beam  volume  con- 
stant.  That  is,  the  beam  volume  must  be  the  same  if  the  variable 
dimension  is  a  constant  or  a  function  of  x.   This  stipulation  will 
become  evident  in  the  derivation  of  the  area  function. 
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DEVELOPMENT  OF  THE  GOVERNING  EQUATION 


With  the  substitutions : 


s=kAG(y-f)  (5) 

(r*  Tat-  <6> 

equations  (1)  and  (2)  become: 

[kifJ]/+s=o  (7) 

f  Ay  -  s'  (8) 

where      =Jj?    and     "  =■  Js 
Solving   equation    (5)    for    \|r : 

l/r=y'-s<r  (9) 

f  =  y  -(S<r)'  do) 

Lfr'^y'-iSG-f  (11) 

Equation    (7)    is : 

£iK+erv+s-  o  d2) 

Substitution  of  equation  (10)  and  (11)  into  (12)  gives  : 

El[^/"Ls"cr-^s/^+EI,Q/,/-5/or]4.s[j-(5IcrO]=0      (#> 

Substituting  equation  (8)  into  (13) : 

El5i'"-crVAi(')-2<r/rA3^El'^"-afA<|]  +  s[i-(EI<r')jr q    (i*) 

Writing  equation  (14)  as: 

A*  +  sB*=o  <14a> 

thus  , 


s--£ 


(15) 


Jl 


where    D:^fO(see   Appendix  B) 
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Differentiating  equation  (15)  with  respect  to  x  gives: 

/>A?B,a+A"'B*-A*B*  =  0  (16> 

This  form  of  the  equation  does  not  seem  to  have  been  previously 
reported.   It  is  clear  that  this  single  equation  is  fundamental  to 
the  study  of  vibrations  of  non-uniform  beams,  including  shear. 

For  the  case  of  a  prismatic  beam  equation  (16)  reduces  to: 

ei  &  -  ifa&e  +  r A^  =  o 

From  this  point  on,  the  analysis  in  this  thesis  will  focus  upon  and 
be  confined  to  the  case  of  simply  supported  beams  for  which  the  end 
conditions  are  : 

(ike  =  fc  w = <yko  =  &  W  =  °         <i?> 

These  conditions  prescribe  zero  deflection  and  zero  moment  at  both 
ends  of  the  beam.   Solutions  to  equation  (16a)  are  : 

Sfak  gn(x>  sm(uVt-KjO     ^  n.9.l,A,3...  (is) 

gn0O=  s/n  AnX  ;    X,=  "^  <19> 

where  : 

are  the  eigen  frequencies  of  vibration  for  the  uniform  beam.   The 

associated  mode  shapes  are  the  corresponding  g  's. 

n 

Equation  (16)  is  completely  general  and  it  is  clear  that  difficulty 
will  attend  attempts  to  find  useful  solutions  which  also  satisfy 
equation  (17).   Therefore,  in  order  to  make  the  equation  more  tract- 
able, a  specialization  is  made  at  this  point.   The  area  function  is 
limited  to  widthwise  variations  so  that  the  moment  of  inertia  is 
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directly  proportional  to  the  cross  sectional  area  and  the  shape  of 
the  cross  section  is  limited  to  those  for  which  the  shear  coefficient 
is  constant. 

After  expansion  and  combination,  equation  (16)  becomes: 

+©k,  \j"0+  ^fe  fwy +eec,y*-©f>wij]  ■*■ 

J^^fey-e^J  4        {20) 

$  fea^ae^!]  -J-2*^'! Bay' 

Equation  (20)  is  the  governing  differential  equation  for  the 
vibrating  Timoshenko  beam,  neglecting  rotary  inertia,  when  the 
shear  coefficient  is  constant  and  the  moment  of  inertia  is  directly 
proportional  to  the  cross  sectional  area. 
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SOLUTION  OF  AREA  FUNCTION 

Because  of  the  high  degree  of  non-linearity  of  equation  (20),  it 
does  not  lend  itself  to  closed  form  solution.   Accordingly,  an 
approximate  method  of  solution  is  employed.   It  is  desirable  to 
obtain  a  solution  which  permits  verification  for  a  prismatic  beam. 
With  these  ends  in  mind,  a  perturbation  expansion  of  the  symmetric 
mode  shapes  is  considered. 

^ssmars  +  €  sin  3ap£  +  &sm  5xgL  (21) 

Physical  arguments  suggest  that  the  substitution  of  equation  (21) 
into  (20)  will  generate  a  symmetric  area  function.   Substitution  of 
equation  (21)  into  (20)  yields  a  differential  equation  with  the 
parameter  g  in  it ,  and  we  will  take  e  to  be  small.   That  is,  a 
perturbation  method  of  solution  will  be  employed.   This  solution 
will  be  a  perturbation  from  the  fundamental  mode  of  the  prismatic 
beam  solution.   If  the  variable  area  function  is  appropriately 
defined,  the  perturbation  method  will  provide  a  link  to  the 
prismatic  beam  solution.   This  link  will  be  effected  by  setting 
the  perturbation  parameter  equal  to  zero  in  the  expression  for 
the  variable  area.   Let: 

A*  Ao  +  eA,  ^e^A*  (22) 

U>*=LQ*  +€uJ?+€*'Ujt  (23> 

In  equation    (22)   A      is    presumed   to   be  constant   since    the   solution 

we   are   seeking    is   a    perturbation  with  respect    to   the    solution   for 
a    prismatic   beam. 


15 


It  is  expected  that  the  solutions  of  equation  (20)  will  be 
valid  only  for  small  values  of  e  since  they  will  be  perturbations 
from  the  prismatic  beam  solutions. 

When  the  perturbation  expansions  (21),  (22)  and  (23)  are  sub- 

2 
stituted  into  equation  (20),  area  terms  of  higher  order  than  e 

are  generated.  These  terms  may  be  neglected  because  the  perturba- 
tion, in  this  case,  considers  only  terms  of  second  or  lower  order. 
For  example,  consider  the  third  term  of  equation  (20): 

A'aA"=  [e  A>63  Al]1  [e  Af+  eJ  A'Q 

-.  [«Vf  *+  &  A* + 2  e3A',  A 'J  [e  A?+  t»Af] 

3 
The  lowest  order  of  e  generated  here  is  e  ,  therefore,  all  terms  are 

2  2 

of  higher  order  than  s   and  the  term  A'   A1   may  be  neglected. 

When  equations  (21),  (22)  and  (23)  are  substituted  into  equation 

(20)  and  common  powers  of  e  are  equated,  three  linear  differential 

equations  in  A- ,  A..  ,  and  A9  are  generated.   The  equation  for  A  is 

simply  the  constant  cross  sectional  area  case. 

E  AX,  g  +  9  Aof  "£y»  -  fAoug  v0  =  O 

For  a  simply  supported  beam: 

£J?a      SCffg 


0)_   is  the  leading  frequency  term  in  the  perturbation  expansion  (23). 

The   first  order  of  e  generates  a  linear  differential  equation 
in  A1(x)  . 

ATtef^V-eK,^]  -A'feer««£*  -e^V 
-ec,**-  efc-c,v,?j+A;[2ec^"^A,  [aec,^ « 

-a^oj^+^e^o1^^  -A0  CEc,^f-f^oa\/, 

This    is   a    linear   differential   equation  with   known  variable   coefficients 


(24) 
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The  only  restriction  on  the  total  area  function  is  that  the  volume 
of  the  beam  must  remain  constant.   Therefore: 

jAAdx  =  VoUme«AoL+efAA,d*  +ea(AAadx  (25) 

The  perturbation  parameter  e  is  arbitrary  thus  defining  a  class  of 
areas.   Therefore,  to  keep  the  beam  volume  constant  for  all  e ,    it 
is  sufficient  to  require  that: 

5AA,dx=SAA,dx=0  C"> 

Assuming  that  A, (x)  is  piecewise  regular  and  symmetric  it  may 
be  represented  by  a  Fourier  series  of  the  form: 

Equation  (26)  restricts  the  arguments  of  the  Fourier  expansion  to 


even  powers    of   TDL  . 


Sm  CTTrQ  I  CaO,  l)l,3r*. 
COS  ITT-  I  =  o  ;   L  =  0,A,*fr.. 

Therefore,  for  \  A^^s  Q  ,  it  is  sufficient  that  the  arguments 
of  the  sine  and  the  cosine  functions  in  the  Fourier  series  be  even. 
Substituting : 

Nf,=  S»K\  Big.  j o4=  Ig. 

into  equation  (11)  gives  : 

Examination  of  equation  (27)  shows  that  only  the  cosine  terms 
in  the  Fourier  expansion  remain  in  the  solution  for  A,(x)  .   if  the 
sine  terms  were  allowed  to  remain  they  would  generate  products  of 
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cos   a  and   cos    not  and    products    of   sin  a  and   sin   na.      These    products 
can   only   be    represented   as    cosine    functions.      The    right-hand   side 
of  equation    (27)    contains    only  sine    terms.      Therefore    the    coefficients 
of    the    sine    terms    in   the   Fourier   expansion  must   all   be   zero.      On   the 
other   hand   the    cosine    terms    in   the   Fourier   expansion  generate    only 
products    of   sine   and   cosine  which  may  be    represented   entirely  as    sine 
functions    through    trigonometric    identities.      Therefore: 


A,M=  ll  MnCOSKoC 


k'  *=o 


Substitution  of  this  expansion  into  equation  (27)  gives: 

oo 

t^Sin(<*-rM)l+(^J|<^^  -     (28) 

There  are  only  two  values  of  n  which  can  satisfy  this  equation: 
n  is  taken  equal  to  +  2  because, 

COS  ZoL  =  cos  (-^°0 

Substituting  n  =  2  into  equation  (28)  gives  : 

M,Z,sir\3oc  +  KZa  smoc=A0jjksS!n  3^  +  k4smoJ| 

Expansion  of  the  terms  K,  and  Z„  show  them  to  be  identically  zero 
(see  Appendix  A).   Therefore,  equating  coefficients  of  common 
arguments  of  the  sine  function  gives: 

M,=  -I?  A0 
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The  unknown  frequency  is  contained  in  the  term  K5. 
Therefore : 

CjS  O  (30) 

Thus,  the  first  eigen  frequency  has  no  contribution  from  the  first 
order  e  solution. 

The  second  order  terms  of  e  in  equation  (20)  produce  a  linear 
differential  equation  in  A~ (x)  with  known  variable  coefficients. 

AaK,cos<x  +  Atl^sincx-Aikscos*.  -AAKHsir\oc= 
f,sm<*  +  fa, sin  3«-*-^3sin5oc  (3i) 

Equation  (31)  is  similar  to  equation  (27).   To  maintain  constant 


volume,  it  has  been  required  that 
L 

A,(x)ck=G 


£ 


o   A> 


The  solution  of  A  (x)  ,  therefore,  will  be  of  the  form: 

Aa60*  J^flncosnoc  <32> 

rv=o 
Substitution  into  the  differential  equation  (31)  gives 

The  value  of  n  =  +  2 ,  +4  are  taken  because  cos  x  =  cos  (-x). 

Aa0<)=  MaCOS«2*+M3COsMcX  (33) 

2 
This  is  the  part  of  the  area  function  whose  coefficient  is  e  . 

Substitution  of  equation  (33)  into  equation  (31)  gives: 

MaZasin<x  +  [KZ,+  M3Z3]sm  3oc  +  IM^Z4  sin  5"oc= 

^,sm*  hh  ^asir\3oc  +  ^3sin  ^04  (34) 

Z9   is  identically  zero  (see  Appendix  A),  therefore: 


M*Za=f,  =  o 
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2 
The  frequency  <jn    is  contained  in  the  term  ^  . 

?,=  fH-Adu£[p+efte)*J=o 

U;* -£i =■  (35) 

To  find  the  unknown  amplitudes,  coefficients  of  sin  5a  and  sin  3a 
are  equated. 

At  this  point  the  class  of  areas  in  terms  of  the  parameter  e  is 

2 
determined  to  the  second  order  of  e,  to,,  is  the  contribution  to 

2 
the  first  eigen  frequency  from  the  e   solution. 

The  frequency  of  the  first  mode  of  vibration  has  also  been  deter- 


(36) 


mined    to   the    second    order    of  e. 

,     gcffrWVA,  (37, 

Equations  (27)  and  (31)  are  third  order  ordinary  differential 
equations.   They  both  have  three  boundary  conditions  imposed  on 
them.   For  equation  (27)  the  boundary  conditions  are: 

x=o    ;    A,=  M, 
x=L    ;   A.- Mi 

L 

A,dx  =  0 


The   conditions    on  A„    are : 

JAadx=G 
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The  solutions  obtained  for  A.,  and  A2  satisfy  both  the  differential 
equation  and  the  boundary  conditions  in  both  cases.   These  solutions 
are  then,  unique  [13], 

The  class  of  areas  obtained  here  are  symmetric  about  the  midpoint 
of  the  beam  for  all  values  of  e.   For  e  =  0  it  may  be  seen  that  the 
solution  for  the  prismatic  beam  is  obtained. 
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SOLUTION   OF   EIGEN  FREQUENCIES   AND  MODE   SHAPES 

Extensive    algebraic   manipulations    have    been   carried    out    through- 
out   the    following   work.      For    the    sake    of   clarity   of   presentation,    many 
long   algebraic   expressions   are    represented   by   symbols.      The   expansion 
of    these    symbols   may  be    found    in  Appendix  A. 

Using    the    results    obtained    for  A,    and   A„  ,    it    is    now   possible    to 
obtain  equation    (20)    in   the    form: 

For   the    simply   supported   beam.,    the    boundary   conditions    are: 

The    procedure    in   solving   equation    (38)   will   be    similar   to   the 
procedure    used    in   obtaining   the   variable   area    functions   A,    and  A„  . 
That    is,    three    separate   differential   equations    in  y         ,    y      ,    and 
y        will   be    obtained    by   equating    like    powers    of   r.    on   either   side 
of   the   equation.      Solutions    to   these   three   equations   will  give   all 
the    eigen    frequencies    and   deflection  modes    in   terms    of   x   and    the 
perturbation   parameter  e. 

The   notation   used   will   be   as    follows: 

^rv*  ^rvo  +e\/ru+ea,}na  (39) 

UOn  =  UJ>no  +  €<^>*>   +  &U)njL  (40) 

The    first    of    the   double    subscripts    refers    to   the   eigen   frequency 
or  deflection  mode.      The    second   subscript    refers    to   the    order   of  e. 
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The  zero  order  of  e  gives  the  constant  cross  section  equation. 

EC,\("+9|»^>2  -f<^>™\j«0  =  O  (41) 

The  mode    shapes    are    of   the    form: 

N^no  =  T3no  Sin   XrvX        ^      "^r\'JQU  (*2) 

The  amplitude,  B  _,  may  be  determined  from  initial  time  conditions. 
Substitution  of  equation  (42)  into  equation  (41)  gives: 

BnoSin  *aX[gC.  l$-&?(*>&  &-f^£j=  O 

Therefore  : 

/     )S  —  fc-      '  ***■  (43 1 

^o-  fD+eXy 

This  result,  for  n  =  1,  checks  with  the  earlier  result  in  Section 

IV. 

Equation  (24)  may  now  be  expanded  into  a  differential  equation 

in  y    .   By  equating  the  coefficient  of  e  in  equation  (38)  to  zero, 
nl 

we  get: 

-  Qt  cos  2<*  Si  a  rv<* 

Applying  trigonometric  identities : 

Q,>.    +QaN/n!  -Q3^n.  =  QqSm  n*  (44) 

where , 

If  Q,  were  non-zero,  its  contribution  to  the  solution  of  equation 
(29)  would  be  non-periodic  [7].   Therefore,  its  value  must  be 
chosen  such  that  it  vanishes  for  all  values  of  n. 
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A  Fourier  series  is  assumed  for  y  , .   The  coefficients  of  all  the 

nl 

cosine  terms  in  this  series  must  be  zero  because  the  right  hand  side 
of  equation  (44)  contains  only  sine  terms.   After  substitution  of  the 
Fourier  sine  series  into  equation  (29),  only  three  sine  terms  have 
non-zero  coefficients. 

Nrvi^  Bru  smCi-aW^Cm  sin  riot  +DnrsinCi+a)oC       (45) 

The   requirement    of  Q.    =0   for   all   n  makes    C    ,    =  0   for   all   n. 

4  nl 

Substitution  of  the  assumed  shape  for  y  n  into  the  differential 

nl 

equation  (44)  gives : 

where  , 

When  n  =  1,  the  sin(n-2)a  term  becomes  -sin  a. 
Therefore : 

Expansion  of  the  quantity  ?|Qs+Qtl  ^or  n  =  1  shows  that  it  is 
identically  zero  (see  Appendix  A).   Therefore: 

f; 

Qm=0 

2 

This  confirms  the  earlier  finding  in  the  A1  (x)  solution  of  a),  =  0. 
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To  maintain  a  periodic  solution,  Q,  must  equal  zero  for  all 
other  values  of  n.   Therefore: 

The  amplitudes  B  ,  and  D  ,  are  found  by  equating  terms  with 
v  nl      nl 

the  same  sine  function  arguments.   For  n  =  1: 

Q,p*-<Sfcf£-Qa 


Ui* 


Qs-Q<, 


The  expansion  of  D,..  shows  it  to  be  identically  equal  to  unity 
This  checks  the  earlier  substitution  of: 

.  n>l 


For   other   n: 


Table    2 


n 

Amplitudes 

Frequency 

KJAI 

Dru 

CiJni 

1 

o 

/ 

O 

n>l 

-CQs+Q*! 

Tq^-qJI 

o 

^^-Q^-Q^i 

afotf-Qktf-flfl 

a)    ,    is    the   contribution  of   the    first    order   e    solution   to   the    overall 
nl 

2 
eigen   frequency  a)    . 
J      n 

B    ,    and   D    ,    are    the   amplitudes    of   the    sine    functions   which 
nl  nl  v 

contribute    to   the    overall   equation   for    the   mode   shapes. 
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where , 

The  second  order  e  equation  in  the  unknown  y    is : 

n2 

+  <S3  5m(n+^)oc  4-  64  Sir\Cr\-^)«  +  4^iaaa  (46) 

The  assumed  shape  of  y  ~  is  therefore: 
The  term  6,  is: 


^5-=  BAofD+QX^Ju;^+4 


As  was  the  case  with  the  first  order  c  solution,  this  term  may  not 

be  allowed  to  persist  if  the  solution  is  to  be  periodic.   This 

condition  defines  the  frequency  and  demands  that  G  „  =  0. 

n2 

For  n  =  1,  the  6~  sin(n-2)a  term  becomes  -  6„  since. 
Therefore  : 


oa-6( 


Similarly  for  n  =  2,  the  5,  sin(n-4)a  term  becomes  -  6,  sin  2a 
and  combines  with  the  &s  sin  2a  term.   Therefore: 

B^pLl  +  HG^^jA'o 

The  amplitudes  B  0 ,  D  -  ,  E  0 ,  and  F  _  are  found,  as  in  the  first 
nZ    n£        nZ       nZ 

order  g  solution,  by  equating  coefficients  of  the  sine  functions 
with  like  arguments.   Care  must  be  taken,  however,  noting  combina- 
tions for  certain  values  of  n.   Substitution  of  the  assumed  form 

of  y   into  the  differential  equation  (46)  gives  : 
n2 
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Bvi[Qlp'1Qap/-Q373i^c^«H)o<=j^,^(a(h4i)<+6ism(n--a)oc 

where , 
For  n  =   1 

[Q,(?F)M-  0*  GST-  Qs][B.a-  R  J  •&-  &»J  i* 

The   expansion   of   6,    -    6,    shows    it    to   be    identically   zero   for    n  =   1 
(see   Appendix  A).      The   expression  K^(^)  "QaljFj  "  Q3J 


is  non-zero, 


Therefore  [B..„  -  F  ~],  which  is  the  coefficient  of  sin  3a  for  n  =  1, 


is  zero. 


2 
From  the  earlier  calculation  of  0),^,  we  know  that  [6c.  -  6~  ]  =0. 

Therefore  : 

Expansion  of   'y  ("Q  q,VL  Q  Q'^— O  I  A   shows  it  to  be  identically 
equal  to  one  for  n  =  1  (see  Appendix  A). 

This  checks  the  earlier  substitution,  the  the  A~  solution,  of: 
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2 
Expansion  of  the  expression  for  go..,,  also  shows  it  to  be  consistent 

2 
for  the  value  go~  found  in  the  A„  solution. 

For  n  -  2  the  sin(n-2)a  term  vanishes.   Similarly,  for  n  =  4  the 

sin(n-4)a  term  vanishes.   Therefore: 

Daa.  =  o 

For  all  other  n. 

c,,_ ^a ,     .  *\^  I 

go  „  is  the  contribution  of  the  second  order  e  solution  to  the  overall 
n2 

■  2 

eigen  frequency  go  . 
n 

B  „,  D  „,  E  „,  and  F  „,  are  the  amplitudes  of  the  sine  functions 
n2    n2    n2       n2  r 

which  contribute  to  the  overall  equation  for  the  mode  shapes, 
where  , 

The  mode  shapes  for  the  class  of  areas  derived  in  Section  IV  are: 
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"Brno  £Ci  X*.    ^g^A-rt 


The   eigen   frequencies   are: 

LJ*.  '  jBaoEC,^^^   ,x„  (48) 

where 

B    is  an  amplitude  which  depends  on  the  initial  time  conditions 
no  r 

of  the  problem. 

By  setting  e  ~  0  in  equations  (47)  and  (48),  the  solution  is 
reduced  to  the  constant  cross  section  solution.   This  is  expected 
because  when  g  =  0: 

A^A0=  constant 

By  setting  n  =  1  in  equations  (47)  and  (48),  the  values  of 
and  <d,  are  : 

These  results  are  identical  with  the  results  and  substitutions  from 
Section  IV,  the  solution  for  the  area  function.   This  checks  the 
great  amount  of  algebra  which  was  performed  in  obtaining  formulas 
for  the  eigen  frequencies  and  mode  shapes. 
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CONCLUSIONS 

A   fundamental   non-linear   ordinary  differential   equation,    equation 
(16),    has    been  derived    for   the    flexural   vibrations    of   a    non-uniform 
beam,    including   shear. 

Equation    (20)    is    the    particular    form   taken  by  equation    (16)    under 
the   specializing   assumptions    that   moment    of    inertia    is    directly   pro- 
portional  to   cross    sectional   area   and    that    the    shear   coefficient    is 
constant . 

Tables    2   and   3   contain  approximate    solutions    of   equation    (20)    for 
all   eigen   frequencies   and  mode   shapes    in   terms    of   a    perturbation 
parameter   e.      The    perturbation  solutions    are   accurate    up   to   the 
second    order   of   e. 
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EXTENSIONS 

The  preceding  results  could  be  applied  to  beam  optimization  design. 
Using  the  formulas  of  the  preceding  sections,  the  deflection  may  be 
obtained  as : 

It  may  be  desirable  to  minimize  the  deflection  y.   For  example, 
minimization  of  the  amplitude  of  vibration  could  cut  down  on  noise 
level.   Bending  stress  is  a  function  of  the  deflection  and  this 
quantity  could  also  be  minimized.   This  minimization  process  could 
be  approached  through  iteration  on  the  perturbation  parameter  e. 
A  value  of  g  would  be  chosen  and  y(x,t,e)  would  be  evaluated  at 
increments  over  the  beam  length  over  appropriate  time  intervals. 
The  maximum  value  obtained  from  this  process  would  be  saved.   This 
process  would  then  be  repeated  for  other  values  of  e.   The  minimum 
value  of  the  maximums  noted  for  each  iteration  would  give  the  optimum 
va  lue  for  £  . 

It  must  be  noted  that  the  equations  developed  are  valid  only  for 
small  values  of  e. 

Another  possibility  for  further  study  would  be  the  application 
of  this  method  to  the  forced  vibration  case  and  subsequent  optimization 
of  design  for  deflection,  stress,  or  frequency. 

An  experimental  investigation  to  compare  with  theoretical  results 
would  be  useful.   Problems  would  probably  arise  in  the  manufacture  of 
of  the  beam  itself  and  the  development  of  devices  sensitive  and  accurate 
enough  to  detect  variations  from  the  solution  for  a  prismatic  beam. 
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APPENDIX  A 


EXPANSION  OF  SYMBOLS 

This  Appendix  contains  the  expansion  of  the  symbols  used  in  the 
main  body  of  the  thesis  and  those  symbols  used  in  Appendix  D.   Symbols 
equal  to  zero  are  shown. 

The  expansion  of  the  symbols  used  in  the  area  solutions  follows. 

Ki=-.aec,GE)3 

fc,=  j  [fecffi*  -  p^  -  ©r  ^^)aJ » o 

Ks--  ecW- r«>t  -Qr^t&f 

Kr.  -  &eojKif  -  r*>* = °    •  c^'a = ° } 

K,o-  3  k*-  o 

k,s=  -3©(=w  -e*pc*->*(?)*= o 

U=  ©  ecm*  -  3e  P  ^a  -  ©WW*-  ^ eft?)* 
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2, 
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The    following    symbols   were    used    in   the    solution    for    the   eigen 
frequencies    and   mode   shapes. 

On^^raZp  +©f  A*]  a;*  -o   |  co^=o  ,  k*=- not 

Qs  -  Bno  W.«H©P«&  JUBtf  -  0EC,  £  6?)*-  2 fc,  XI] 
QlM  -  Bru  [jETC,  py-  f  <m£>  -  ©f  ^ifo  p4I 

Qn-Ba.^ECp3] 
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Qas  = 
Qa-1 

Q« 

4>.= 

$*- 

4 


=  Bm  L-  ©  EC,  pS+  ©^  C6$  pi 
-  Dm  [r  ©£C,  Q2V  0~f  </j*\  <jj[ 

=  B*o  |©  you; Fi  AkI  -  O 

-  AoMaftS^Qql 

-  A0  Kb  [M1  Q^l 

AoM3[C?)Q?+ffl3Q.o] 
AoQ.. 


A*  Q.3  ,       _ 

A2  W,  fcQiM  -TO*QaHJ 

AJN.CaQw-aB^aJ 

AJ/V.lXW^aj-^JQn] 

AjN^^^Q^+m1^^ 


'N 


^  A:w, 


CWo»-(WQ«. 
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6,=  iD*,+<J>3  +  4>-r] 
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APPENDIX  B 


THE  CASE  WHEN  B*  IS  ZERO 


On  page  10  an  assumption  that  B*  was  not  equal  to  zero  was  made 
Suppose  B*  did  vanish: 

B*--|-[EI<r']  =  0 
[Elcr'J'-j 

Integrating  with  respect  to  x: 


If  k  =  constant  and  I  =  C..A 


KG, 


A»exp[-^^Vax^bj] 


This  is  not  a  symmetric  function  about  the  beam  midpoint.   Since  we 
are  restricting  consideration  to  situations  having  symmetry  about 
x  =  L/2 ,  this  case  is  excluded. 
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APPENDIX  C 


NUMERICAL  EXAMPLE 


Steel  Beam 


L=l00m- 


E=3C>x06ps( 

Gi=iax/o^psi 

r.-wow/ft3 


Brio*  lii\-  ,  r\>/\ 


A  computer  program  was  prepared  that  computes  the  first  ten 

eigen  frequencies  and  mode  shapes  for  g  =  .05.   The  program  also 

2 

computes  and  plots  the  variable  area  function  eA,  +  e  A?  for 

€  =  .01,  .02,  .03,  .04  and  .05.   The  first  ten  eigen  functions  for 

e  =  .05  are  also  plotted. 

2 
The  physical  significance  of  eA1  +  e  A«  in  this  case  would  be: 

A-  book  =   A0  +  €A,+  e*Ax 

Therefore,    the  width  would   be   a    function   of  x  while    the    height, 
h,    remains    constant. 
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The    first    ten  eigen   frequencies   are: 


n 

>* 

2. 

LOnx 

^4*£&€u&. 

i 

.332xl04 

0.0 

-.261xl06 

51.6   cycles/sec. 

2 

.530xl05 

0.0 

-.310xl07 

212.7 

3 

.268xl06 

0.0 

.909xl07 

539.2 

4 

.845xl06 

0.0 

.190xl08 

944.9 

5 

.206xl07 

0.0 

.322xl08 

1,462.4 

6 

,425xl07 

0.0 

.488xl08 

2,091.7        " 

7 

.785xl07 

0.0 

.693xl08 

2,832.4 

8 

.133xl08 

0.0 

.941xl08 

3,683.1        " 

9 

.212xl08 

0.0 

.124xl09 

4,642.4 

10 

.322xl08 

0.0 

.158xL09 

5,708.9 

40 


< 

n 
f 


-a 


G-.oi 

ao 


40 


so 


80 


X  0«O 


+ 
< 


-j 


€=o3 

ao 


40 


40 


80 


/co 


XOM 


< 

-3 


€.=.<3S 

^ 

40  60 

XOM 


SO 


/OO 


41 


(in) 


("irs"^    mode  sliape, 


20 


40  6O 


so 


100 


•Decorvd     rv\odz,    sUpe 


100 


j(ii)0 


-J 


I    h  ire)     mode    S  ka  P1 


«-  ■' 


ao 


"°      XOM  6o 


SO 


100 


42 


yin)  o 


fourtk     mode    skape 


IOO 


/(in-) 


f"  i  tl  n    mode  shape 


too 


\j(<*.)o 


O  I  x"t  k    >wode   sKcxpe 


IOO 
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APPENDIX  D 


INTERMEDIATE  FORMS  OF  EQUATIONS 


H$ 


Referring  to  page  25  the  differential  equation  for  A~(x)  is: 

EC,]*-  frffr  -  f  MJy,  te^Y  4  GpusitfJ  +AoA," 

^ec^'l+A.Ar^fc-c^+Ar^-je^^+eecv,;^ 
A.'Afl*?©^?^  +  g©  6C,  y 1 4  Af[-«:r wty-ee:,^? 

-3©f  "J^,  -  3G(=-6o,i\(0+  ©*r<A^y+©*j»tw,a  \/0" +€C,y| 

Substitution  of  A, (x) ,  y  ,  y  ,  and  y_ ,  and  subsequent  reduction 
gives : 

Ao[A^K,o)SPc+A2kasmo(-ALK3Cos^s^is'A^+^as'n3K 

^€»  s  I  ft  ^  cos  3©<  +  fi-7  coscx  coSc5o<  s;n  ^oi-f 
&^5!MCoS  -204-fAcj  3(fto<  SIA"*^o<   -  S/6  co5<^<  3/n  3*<* 

After  further  reduction,  using  trigonometric  identities,  equation 
(31)  is  produced. 
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The  basic  differential  equation  for  y  0  is  similar  to  that  for 

nl 

A?(x).   The  only  difference  is  that  y   is  treated  as  the  unknown. 

*  n^ 

After  substitution  of  A^  A2 ,  ynQS  and  y    the  equation  takes  the 
form : 

Ao[Q»\/S  tQa^n*  -Q3^iJ=  cf),co5<2^  sin  not  +<|>acos4o<  sin  na 
+c|>3s,n  2<*  cos  n<*  -  c|>4  s<n  ^  cosnc*  +45  sum*  f^smpx 
+^7sm^x  -4%cos^°<s,r>P*  "  4«i  cos  ^  smcjx 
+cj),o5fnAx  cospx  +^,i  smPvoi  coscjx  +^acos^sm^  cosa* 
-KJ>i3  co5^o(  sm  not  +cj>,w  sin^A^sinnoc 

Application  of  trigonometric  identities  to  this  equation  produces 
equation  (46). 
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